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Abstract — This paper presents Chemo-tactic PSO-DE 
(CPSO-DE) optimization algorithm combined with 
Lagrange Relaxation method (LR) for solving Unit 
Commitment (UC) problem. The proposed approach 
employs Chemo-tactic PSO-DE algorithm for optimal 
settings of Lagrange multipliers. It provides high-quality 
performance and reaches global solution and is a hybrid 
heuristic algorithm based on Bacterial Foraging 
Optimization (BFO), Particle Swarm Optimization (PSO) 
and Differential Evolution (DE). The feasibility of the 
proposed method is demonstrated for 10-unit, 20-unit, 
and 40-unit systems respectively. The test results are 
compared with those obtained by Lagrangian relaxation 
(LR), genetic algorithm (GA), evolutionary programming 
(EP), and genetic algorithm based on unit characteristic 
classification (GAUC), enhanced adaptive Lagrangian 
relaxation (ELR), integer-coded genetic algorithm 
(ICGA) and hybrid particle swarm optimization (HPSO) 
in terms of solution quality. Simulation results show that 
the proposed method can provide a better solution. 

Index Terms — Lagrangian Relaxation, Particle Swarm 
Optimization, Differential Evolution, Bacterial Foraging 
Optimization, Unit Commitment 

I. Introduction 

Unit commitment (UC) is used to commit the 
generators such that the total production cost over the 
predicted load demand for scheduled time horizon is 
minimized considering the spinning reserve and 
operational constraints of generator units [1], [2]. Unit 
commitment is a high dimensional, nonlinear, non- 
convex, mixed-integer combinatorial optimization 
problem. Priority list method, integer programming, 
dynamic programming (DP), branch- and-bound 
methods, mixed-integer programming, and Lagrange 
relaxation (LR) are a few methods developed up to now 
for solving UC problem. 

In priority list method, priority of units is determined 
from full load average production cost of the unit. The 
method is simple but the quality of solution is low. 
Priority list of units is also considered in Dynamic 
programming method. In spite of many advantages 
such as ability to maintain solution feasibility, for 



higher dimension problem the problem size increases 
rapidly with the number of generators and this requires 
enormous computation time and large memory space. 
Branch-and-bound and mixed integer linear 
programming method also requires large computation 
time and memory space [3], [4]. Lagrange relaxation 
for UC problem is advanced than dynamic 
programming due to its faster computational time. The 
solution to Lagrange Relaxation for UC problem 
depends on the updating of Lagrange multipliers; hence 
it suffers from solution quality problem. This paper 
proposes a new hybrid heuristic method for solving UC 
problem. The proposed method is developed in such 
way that a Chemo-tactic PSO-DE optimization 
technique is applied to update Lagrange multipliers and 
this improves the performance of LR method. To 
illustrate the effectiveness of the proposed method, it is 
tested and compared to the conventional LR [5], GA 
[5], EP [6], GAUC [7], ELR [8], ICGA [9] and HPSO 
[10] for 10-unit, 20-unit, and 40-unit respectively. 

In 2001, Prof. K. M. Passino proposed an 
optimization technique known as Bacterial Foraging 
Optimization Algorithm (BFOA) based on the foraging 
strategies of the E. Coli bacterium cells [11]. Until date 
there have been a few successful applications of the 
said algorithm in optimal control engineering, 
harmonic estimation [12], transmission loss reduction 
[13], machine learning [14] and so on. Experimentation 
with several benchmark functions reveal that BFOA 
possesses a poor convergence behavior over multi- 
modal and rough fitness landscapes as compared to 
other naturally inspired optimization techniques like the 
Genetic Algorithm (GA), Particle Swarm Optimization 
(PSO) and Differential Evolution (DE). Its performance 
is also heavily affected with the growth of search space 
dimensionality. In 2007, Kim et al. proposed a hybrid 
approach involving GA and BFOA for function 
optimization [15]. The proposed algorithm 
outperformed both GA and BFOA over a few 
numerical benchmarks and a practical PID tuner design 
problem [16-19]. 
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II. NOMENCLATURE 

Generation output power of unit i at hour t. 



Status of unit at hour t (on=l, off=0). 



p; 

ST t : Startup cost of generator i at hour t. 

F i {P j ) : Generator fuel cost function in quadratic 
form 

F l {p;) = a i (p;f+b l p;+c l ($/h) 

F{P i ,U t ) : Total production cost. 

N : The number of generator units. 

N s : Swimming length Ns is the maximum 

number of steps taken by each bacterium 
when it moves from low nutrient area to 
high nutrient area. 
The number of bacteria in the population. 

The number of chemo-tactic steps. 

The number of hours. 

The best previous position of k 

particles. 

The index of the best particle among all 

particles in the group. 

The scale factor 

The Cross over probability 

Acceleration constants 
rand, Rand ; Random value in the range [0, 1] 
^load '• Load demand at hour t. 

R* : Spinning reserve at hour t. 

-S.min : Minimum generation limit of generator i 

^•,max : Maximum generation limit of generator i 

Minimum up time of generator i 

Minimum down time of generator i 

The unit's hot startup cost. 

The unit's cold startup cost. 

The cold start hour. 

Lagrange multipliers 

Velocity vector of 2! in t hour for k 

bacterial population. 
Velocity vector of [it in t hour for k 

bacterial population. 

III. PROBLEM FORMULATION 

From the definition of UC mentioned above, the 
objective of UC problem is to minimize the production 
cost over the schedule time horizon (e.g., 24h). The 



S 

N c 

T 

pbest 

gbest 

F 
CR 

c c 



i,up 

T 

i,down 

HST. 

CST t 
T 

i,cold 

VX(k) 
Vju'(k) 



objective function can be formulated mathematically as 
an optimization problem as follows: 

T N 

f(p; ,uj) = ££[^ {p; ) + st; (i - ur w; (i) 

subject to 

(a) Power balance constraint 



PL-t p ' u '=° 



!=1 

(b) Spinning reserve constraint 

jV 

PL+* t -Y t p i*~ ut t*° 

!=1 

(c) Generation limit constraints 

P. ■ U* < P' < P u' 

/,min / — i — /,max ; 

(d) Minimum up and down time constraints 

1 V ■* i,on "^ •* i,up ' 

U' = \ if T „ <T., , 

i J i, off i, down ' 

or 1, otherwise, 

(e) Startup cost 

, _ \HST t , if T idown < T Uoff < T jcold +T idown 
ST, — 



[CSf, ifT ii0ff >T UcM +T Umm 



(2) 

(3) 
(4) 

(5) 
(6) 



IV. OVERVIEW OF PARTICLE SWARM 
OPTIMIZATION 

The Particle Swarm Optimization was introduced by 
James Kennedy and Russell C. Eberhart in 1995 [20]. It 
is similar to other evolutionary computational 
techniques those are GA and EP in that PSO initializes 
a population of individuals randomly. These 
individuals are known as particles and have positions 
and velocities. In addition, it searches for the optimum 
by updating generations, and population evolution is 
based on the previous generations. PSO is motivated 
from the simulation of the behavior of social systems 
such as fish schooling and birds flocking. In PSO, each 
particle flies through the problem space by following 
the current optimal particles. Each individual adjusts its 
flying according to its own flying experience and its 
neighbors flying experience. Based on its own thinking 
the particle attracts to its best position ( pbest ) the 

particle changes its velocity based on the social- 
psychological adaptation of knowledge that is the 
particle attracts its previous best position among the 
group ( gbest ), here so the velocity of the particle is 

updated to new position according to its modified 
velocity using the information. 

a) The current velocity 

b) The distance between the current position and 
pbest . 

c) The distance between the current position and 
gbest . 



51 



©2010 ACEEE 
DOI:01.ijepe.01.02.10 



4eACEEE 



ACEEE International Journal on Electrical and Power Engineering, Vol. 1, No. 2, July 2010 



V v (iter + 1) = w * V tj (iter) 

+ C, *rand l *(pbest i/ — x t (iter)) (j) 

+ C 2 * rand 2 * (gbest i — x tJ (iter)) 

where 

W = inertia weight factor 

Xy (iter) - current position of i particle in j 

dimension at iteration (iter) 

Vyiiter) — velocity of /'* particle in j' dimension 

at iteration (iter) 

Appropriate selection of inertia weight provides a 
balance between global and local explorations. The 
updated position of each particle is expressed as 

Xy (iter + 1) = Xy (iter) + V tj (iter + 1) (8) 

The new positions are calculated repeatedly till a pre- 
specified maximum number of iterations are reached. 

V. OVERVIEW OF DIFFERENTIAL 
EVOLUTION OPTIMIZATION 

Differential Evolution (DE) is a new floating point 
encoded evolutionary algorithm for global optimization 
developed by Price and Storn in 1995 [21]. It is similar 
to other evolutionary computation technique and also 
starts with a population of NP. It is a D-dimensional 
search variable vector. A special kind of differential 
operator is used to create new offspring from parent 
chromosomes instead of classical crossover or 
mutation. In each generation to change each population 

member Xi, a Donor vector V. is created. Three 

parameter vectors r\, rl, and rl are chosen in a random 
fashion from the current population. A scalar number F 
scales the difference of any two of the three vectors and 
the scaled difference is added to the third one. Donor 

vector for j component of k population at iter 

generation is expressed as 
v kj (iter + l)=X rlJ (iter) 

+ F *(X r2 f (iter)- X r3J (iter)) {9) 

CR is called "Crossover" constant and it appears as a 
control parameter of DE. The crossover is performed 
on each of the D variables whenever a randomly picked 
number between and 1 is within the CR value. The 

trial vector U k . is outlined as 

u h . (iter) = v k . (iter) if rand(0,l) < CR 

u k . (iter) = X k . (iter) otherwise 

VI. Overview of Chemotactic PSO-DE 

PSO and DE are excellent heuristics like other 
evolutionary algorithms. Practical experiences suggest 
that they reach stagnation after certain number of 
generations as the population is not converged locally, 
so they will stop proceeding towards global optimal 



solutions. Chemo-tactic process of bacterial foraging 
describes the movement of an E.coli cell through 
swimming and tumbling via flagella. It can move in 
two different ways, it can swim for a definite period of 
time or it can tumble and it alternates between these 
two modes of operation for the entire lifetime. For 

example if k denotes the fc th bacterium and Ck is the 

size of the step taken in the random direction specified 
by the tumble. Then the movement of bacterium is 
represented by 



k =O k +C k * Delta k 



Delta, 



A, 



(11) 



where A is a unit length vector in the random direction 

In this article we come up with a hybrid optimization 
technique, which synergistically couples the BFOA 
with the PSO, DE and Lagrange multipliers method. 
The proposed algorithm performs local search through 
the chemo-tactic movement operation of BFOA 
whereas the global search over the entire search space 
is accomplished by a PSO-DE operator. In this way it 
balances between exploration and exploitation enjoying 
best of both the worlds. 

VII. METHODOLOGY 

Since UC was introduced, several methods have been 
used to solve this problem. Among those methods, LR 
seems to be the most appropriate one [5], [8]. This 
method solves the UC problem by relaxing or 
temporarily ignoring the coupling constraints, power 
balance and spinning reserve requirements, then 
solving the problem through a dual optimization 
procedure. The dual procedure attempts to reach the 
constrained optimum by maximizing the Lagrange 
function with respect to the Lagrange multipliers. 
L(P, U , A, M) = F(P' , U ' ) 



+jyt\n ad -jy;u; 



(12) 



-I> 



-R' 



^m^ i, max ( 



while minimize with respect to the power output and 
generating unit status, that is 
generating unit status, that is 



q*(A,jU) 



- max q(A., ju) 



13) 



where 

q(A,fj)=m\nL(P,U,A,fJ) (I 4 ) 

p'i fJ] 

Substitute (1) in (9), the Lagrangian function becomes 
\Fi(P') 

l = > '> ' j +5*7;.' (i — u;- ] )] U) 

-XP'U] -jUP im ^U' t I (15) 



+Yi^PL d +U(Kaa +*')) 
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Equations (2) and (3) are the coupling constraints 

across the units. The coupling constraints in the above 

equation are temporarily ignored and then the 

minimum of Lagrangian function is solved for each 

generating unit, without regard for what is happening 

on the other generating units. 

* t \[F : (P') + ST'(l-U'- l )]U'] 

mm L(P,U, I, n) = Y minYT (16) 

Therefore the minimum for each generating unit over 
all time periods is: 

mmqM = ^ mm ^-xm-^u; } (17) 

subject to constraints (4) and (5). This problem is 
solved through a two-state dynamic programming 
problem in two variables for each unit. On the other 
hand, in order to maximize the Lagrange function with 
respect to the Lagrange multipliers, the adjustment of 
Lagrange multipliers must be done carefully. Most of 
research works use sub-gradient method to achieve this 
task. In this paper, we use the Chemo-tactic PSO-DE 
optimization to adjust the Lagrange multipliers and 
improve the performance of Lagrange relaxation 
method. 
This is solved as a dynamic programming problem in 

one variable with two possible unit states ( U t — or 
1). 




t=0 t= t=2 1=3 

Fig'. 1 Unit on/off status 

The dynamic programming part must take into account 
all the start up costs as well as the minimum up and 
down time for the generating units. 

At U t =0 state the value of function to minimize is 

zero. At U j =1 state the function to be minimized is 

min[i^. (P' ) - XP. ] . The minimum of this function 
is found by taking the first derivative 

-^— [ F t i p t ) —-HP,' ] =0 



(18) 



■)t,opt 



The dual power ' is obtained from (19) 
X -b, 



pt.opt __ 



2c ; 



(19) 



ifP!-" P '<P^ then PI = P imm 
if P ;** > P t ^, ihenP;=P t ^ (20) 

if P < P''" p ' < P , then P' = P*'" p * 

J z.min — i — /.max ' t,lt " t ; ; 

The dual power obtained from (20) is substituted in 
on/off decision criterion DC 

DC] =[F i {P<) + ST{(\-U<- i )-XP; -M'P^l (21) 



If DC i < the unit will be committed if it does not 

violate minimum downtime constraint, otherwise the 
unit is decommitted if it does not violate the minimum 
uptime constraint. A step-by-step procedure for the 
proposed method is outlined as follows. 

Section -A: 

Stepi: Compute X and JLl' by Chemo-tactic PSO-DE. 

Step2: Set the unit status by calculating dual power and 
obtain on/off decision DC (21). 

Step3: Calculate the spinning reserve for generated 'U' 
status at each hour of the population. If excessive 
spinning reserve is observed decommit the unit with 
high full load average production cost if it does not 
violate the minimum uptime constraint. Otherwise if 
the unit violates, the unit with next high full load 
average production cost is decommitted. Now set the 
spinning reserve according to the new status and repeat 
until the spinning reserve satisfies. 
Step4: If less spinning reserve than required is observed 
commit the unit with less full load average cost if it 
does not violate the minimum downtime constraint. 
Otherwise if it violates, the unit with next less full load 
average production cost is committed. Now set the 
spinning reserve according to the new status and repeat 
until the spinning reserve satisfies. 
Steps: Economic dispatch for the generated U-status is 
carried out by Lagrangian multiplier method. 

Section-B: 

Stepi: Randomly generate X > [1 > VX > an d VjU with 
in the range for initialized bacterium. Generate unit- 
status U as in section-A and calculate pbest for U, X 
, U 1 , cost and fit (J and JF) and gbest for U, X , H l , J 
and JF. 

Step2: Starting of the chemo-tactic loop (ir=it+l) 
Step3: For each bacteria (k=k+l) 
JFlast =JF(k) 

Step4: Calculate VX , VfX and check for maximum 

and minimum velocity limits. 

VX(k)=VX(ls) 

+ C, * Rand * {pbestZ (k) —Z (k)) (22) 
+ C 2 * Rand * (gbestZ —Z (k)) 

VjU (k) =VjU (k) 

+C, * Rand * (pbestjU (k) —/J (k)) (23) 
+C 2 * Rand * (gbestjU —/J (£)) 

Steps: Move to new position 

SL(k) =JL(k) +C^ *vA.{k) 

jU(.k)=jU{k)^-C M *V/J{k) ( 24 > 

where C x and C ^ are chemo-tactic step size 
Step6: Handle limits 

• ^ / ^ (25) 
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Step7: Generate U-status as in section-A 

Steps: Compute cost and fit - J and JF 

Step9: Swimming 

9a) swim-counr=0 

9b) while (swim-count<Ns) and JF(k) > JFlast 

i. swim-count = swim-count + 1 
ii. JFlast =JF{k) 

X(k)=X(k)+C x *VX(k) 

'"■ A /(k)= M '(k)+C fi *V A /(k) 

iv. Handle limits 
SteplO: Generate unit-status U as in section-A 
Stepli: Compute cost and fit - J and JF 
Stepi2: Differential Evolution 
While(rand<CR) 

Xdef (k) = A'(k) + F* (A' (r 2 ) - X (r 3 )) 

fldef'(k) = fl'(k)+FHfl'(r 2 )-fl'(r 3 )) (26) 
i. Handle limits 
Stepl3: Generate unit status —**> fat as in section-A 

Stepi4: Compute cost and fit - J de f , JFdef 
Stepl5: 

if(JF def >JF)then 

JF = JF deP J = J*f,U = U def ,X = X def ,\i = n def 

Stepl6: Update best location pbest and gbest for JF, 

J, U, X, ft 
Stepi7: if k<S, go to step3 otherwise perform stepl8. 

Stepis: If iter < N c then go to step2 otherwise stop. 

VIII. RESULTS AND DISCUSSIONS 

The 10-unit system data and load demands are taken 
from [5]. The 20 and 40 unit's data are obtained by 
duplicating the ten unit case, and the load demands are 
adjusted in proportion to the system size. In the 
simulation, the reserve is assumed to be 10% of the 
load demand. The control parameters chosen for these 

units are CR=0.7, F=0.1 or 0.2, Ci=0.5, C 2 =0.5, C x and 

CjU ^random number between 0. 1 and 1 . The unit 

status obtained through Chemo-tactic PSO-DE for 40 
unit system is shown in Table I. From simulation it is 
observed that setting up of maximum and minimum 
bounds for X and (t, such that Amax = 30 , 
/Imin = 12 , //max = 15 and //min = 2 , leads 
the system to high-quality convergence. The proposed 
method provides best production cost when compared 
to literature as in Table II with a reasonable 
computation time per iteration of 0.72s for 10-unit 
system, 2.1s for 20-unit system and 5.7s for 40-unit 
system. The convergence characteristics for the 
proposed method (40-unit system) are as shown in 
figure2. 



TABLE I 

Unit status for 40-unit system 



Peri 


40-Unit Status 


od 


1-40 


1 


1 1 1 1 1 1 1 100000000000000000000000000000000 


2 


1 1 1 1 1 1 1 100000000000000000000000000000000 


3 


1 1 1 1 1 1 1 1 00000000000 1 00000000000000000000 


4 


1 1 1 1 1 1 1 10000000101 1 100000000000000000000 


5 


1 1 1 1 1 1 1 10000001 1 1 1 1 100000000000000000000 


6 


1 1 1 1 1 1 1 1 000 1 1 1 1 1 1 1 1 1 00000000000000000000 


7 


1 1 1 1 1 1 1 1 00 1 1 1 1 1 1 1 1 1 1 00000000000000000000 


8 


111111111111111111 1 100000000000000000000 


9 


111111111111111111111111 1000000000000000 


10 


11111111111111111111111111111111 00000000 


11 


1111111111111111111111111111111111 1 10000 


12 


1111111111111111111111111111111111111111 


13 


111111111111111111111111111111 1 100000000 


14 


11111111111111111111111 1 1000000000000000 


15 


111111111111111111 1 100000000000000000000 


16 


111111111111111111 1 100000000000000000000 


17 


111111111111111111 1 100000000000000000000 


18 


111111111111111111 1 100000000000000000000 


19 


111111111111111111 1 100000000000000000000 


20 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 10001 1 1 1 1 1 1 10000 


21 


11111111111111111111111 1 1000000000000000 


22 


111111111 10000001 1 1 1 1 1 1 1 1000000000000000 


23 


11111111 0000000000 1 1 00000000000000000000 


24 


1 1 1 1 1 1 1 100000000000000000000000000000000 



TABLE II 
Comparison OF COST 



Total Cost (S) 


No. of Units 


10 


20 


40 


LR[5] 


565825 


1130660 


2258503 


GA[5] 


565825 


1126243 


2251911 


EP[6] 


564551 


1125494 


2249093 


GAUC [7] 


563977 


1125516 


2249715 


ELR [8] 


563977 


1123297 


2244237 


ICGA [9] 


566404 


1127244 


2254123 


HPSO [10] 


563942 


— 


— 


Proposed method 
CPSO-DE 


563977 


1123297 


2243363 



° 2350000 - 



O 2300000 - 

r- 




1 131 261 391 521 651 781 911 
Iterations 



Fig. 2. Cost convergence characteristics for 40-unit system 
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DC. conclusion 

Unit commitment problem is solved with a new 
methodology Chemo-tactic PSO-DE Optimization 
algorithm. PSO and DE are incorporated in chemo- 
tactic algorithm to update Lagrange multipliers so that 
they are suitable for high dimensional combinatorial 
optimization problem. Results show that as the 
dimension of the problem increases it generates good 
unit status with lower production cost. Hence it can be 
concluded that the proposed method provide lower 
production cost than those of LR, GA, EP,GAUC and 
ELR methods for 40-unit system and also the problem 
can be extended for 60-unit, 80-unit and 100-units as in 
this method better costs are obtained as the dimension 
(number of generators) of the problem increases. 
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